Measuring Z' couplings at the LHC 

Frank Petriello and Seth Quackenbush 

Department of Physics, University of Wisconsin, Madison, WI 53706, USA 



Abstract 



We study the properties of potential new Z' gauge bosons produced through the 
Drell-Yan mechanism at the LHC. Our analysis is performed using a fully differen- 
tial next-to-leading order QCD calculation with spin correlations, interference effects, 
and experimental acceptances included. We examine the distinguishability of different 
models and the feasibility of extracting general coupling information with statistical, 
residual scale, and current parton distribution function error estimates included. We 
extend a previous parametrization of Z' couplings to include parity-violating coupling 
combinations, and introduce a convenient technique for simulating new gauge bosons 
on-peak using the concept of basis models. We illustrate our procedure using sev- 
eral example Z 1 models. We find that one can extract reliably four combinations of 
generation-independent quark and lepton couplings in our analysis. For a Z' mass of 
1.5 TeV, one can determine coupling information very well assuming 100 flT 1 of in- 
tegrated luminosity, and a precise measurement becomes possible with 1 ab _1 at the 
SLHC. For a 3 TeV mass, a reasonable determination requires the SLHC. 



1 Introduction 



Z 1 gauge bosons that appear in U(l) gauge extensions of the Standard Model (SM) are 
the most ubiquitous particles in models of new physics. They appear in grand unified 
theories such as S'O(IO) [1] and E(Q) [2], in Little Higgs models [3], and in theories with 
extra space-time dimensions [4]. They often appear as messengers which connect the SM to 
hidden sectors, such as in some models of supersymmetry breaking [5] and in Hidden Valley 
models [6]. Z' states that decay to lepton pairs have a simple, clean experimental signature 
and can easily be searched for at colliders. Current direct search limits from the Tevatron 
require the Z' mass to be greater than about 900 GeV when its couplings to SM fermions 
are identical to those of the Z boson [7]. 

Since the experimental signature is clean and the QCD uncertainties for inclusive quanti- 
ties such as the total cross section and spectrum have been studied and found to be fairly 
small [8], it is likely that the couplings of a discovered Z' can be studied with reasonable 
accuracy to probe the high scale theory that gave rise to it. Many studies of how to measure 
Z' properties and couplings to SM particles have been performed [9]. In particular, a recent 
study focusing on Tevatron physics introduced a parametrization of the parity symmetric 
couplings of the Z' to SM fermions that allows for a simple comparison between experimental 
measurements and theoretical models [10]. 

We attempt to extend previous studies of Z' coupling extractions in several ways in this 
paper. 

• We perform a next-to-leading order (NLO) QCD calculation of the fully differential 
PP (7? Z, Z') — > l + l~X cross section with all interference effects and spin correlations 
included. We also include realistic LHC acceptance cuts. 

• We study the effect of statistical, parton distribution function (PDF), and residual scale 
errors on important Z' observables such as the total cross section, forward-backward 
asymmetries on and off-peak, and central to forward rapidity ratio. 

• We extend the parametrization of Z' couplings in [10] to include parity- violating cou- 
pling combinations, which can be accessed when differential measurements are made. 

• We introduce the use of basis models when simulating Z' states. These arise from 
the observation that the differential cross section can be written as a product of Z' 
couplings multiplied by functions that depend significantly on the specific Z' under 
consideration only through its mass. These functions depend on the PDFs, matrix ele- 
ments, acceptance cuts, and details of the experimental analysis, but need be computed 
only once for a given Z' mass. They can be obtained by running Z' simulation codes 
for basis vectors in coupling space. This facilitates the separation of the specifics of a 
given Z ! model from the details of QCD and the experimental analysis. The functions 
we introduce encoding these details are extensions of the w u ^ introduced in [10]. 

• We study the extraction of the parity symmetric and parity violating Z 1 couplings at 
both the LHC and SLHC, and quantify the effect of statistical, PDF, and residual 
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scale errors on the accuracy of their determination. We note interesting correlations 
between these errors. Using several example Z' models, we examine how well the LHC 
and SLHC can distinguish between different models. 

We use four example Z' models to illustrate our techniques: three models arising from an 
E(6) unified theory and one from a left-right symmetric model. We find that residual QCD 
scale uncertainties have a negligible effect on the measurement of Z' couplings. Statistical 
and current PDF errors are larger and have an approximately equivalent effect at the LHC. 
PDFs should become more accurately known improve when LHC data is used to constrain 
them. We find that the parity symmetric couplings of a 1.5 TeV Z' should be measured with 
good precision at the LHC. Although some information can be obtained about the parity 
violating couplings at the LHC, a more accurate determination requires super-LHC (SLHC) 
statistics. Measurements of the couplings of a 3 TeV Z 1 require SLHC statistics. 

Our paper is organized as follows. In Section [2] we review the theoretical assumptions 
underlying our analysis and present details of the four example Z' models we consider. In 
Section [3] we present the Z' observables used in our analysis. We discuss the details of 
our calculation and present results showing the effect of statistical, PDF, and residual scale 
uncertainties on the basic Z' observables in Section |4} The bulk of our study is discussed in 
Section |5j We introduce our parity symmetric and parity violating Z 1 coupling combinations, 
and show how to extract them from on-peak measurements using the concept of basis models. 
We also study how well the LHC and SLHC can determine Z 1 couplings using our four 
example models for illustration. In Section [6] we present our conclusions. 

2 Theoretical framework 

We first describe the types of Z 1 models we consider in our study. We assume that the Z' 
couplings are generation independent to avoid large flavor changing neutral currents that 
would restrict Mz* to be 100 TeV or more. We also assume that the generator of the U(l) 
group giving rise to the Z' state commutes with the Standard Model SU{2)l generators. 
This implies that the couplings of the Z' to ul and g^, the members of an SU{2)l quark 
doublet q^, are the same, as are the couplings of the Z' to the members of a lepton doublet 
II. These restrictions leave us with the following five parameters: the coupling to q^, the 
coupling to II, and the three couplings to the SU{2) L singlet states ur, d R , and e#. We 
have absorbed the overall gauge coupling into these five quantities. We neglect possible 
Z — Z' mixing; LEP Z-pole measurements restrict the mass mixing angle to be less than 
approximately 10~ 3 radians [11]. 

We utilize four models representative of Z' models discussed in the literature as examples 
to illustrate the extraction of Z' couplings from LHC data. In particular, we examine three 
possible U(l) Z' bosons originating from the exceptional group E$, and one coming from a 
left-right symmetric model, which can arise from an 5*0(10) GUT. We describe these models 
below and list the couplings of the SM fermions to the Z'. 
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• Eq. E 6 models are described by the breaking chain 

E 6 -> SO(10) x -> 5C/(5) x C/(l) x x C/(l)^ -^SMx [/(l)^ (1) 

where 

Z' = Z;cos/3 + Z;sin/3 (2) 

is the lightest new boson arising from this breaking. In this paper we examine the x 
model (P = 0), the tp model ((3 = 7r/2), and the rj model ((3 = arctan(— a/5/3)). 

• Left-right symmetric models: We also consider a left-right model coming from the 
symmetry group SU(2) R x SU(2) L x U(1)b-l- Left-right models can arise from the 
following breaking of SO (10): 

50(10) -> SC/(3) x 5C/(2)i x 5^(2)^ x C/(1) B _ L . (3) 

The Z' in left-right models couples to the current 

j lr = a LRJm- ^^Jb-l, ( 4 ) 

with ctLR = V( c w9r/ s w9l) ~ 1> anc ^ 9l = e/ sin Ow- In the symmetric case $l = 9r 
and a LR ~ 1.59, using the on-shell definition of sin 2 9 W . The overall coupling strength 
is e/ cos#w We use the symmetric model in our analysis. 

The fermion coupling assignments in these models are summarized in the following table. 
For convenience, we factor out an overall e/ cosQw- We assume that the couplings for E 6 
models retain their GUT-scale relations to the EM coupling down to the Z' scale to good 
approximation. Assuming a different value for the overall coupling may have some impact 
on the total cross section. However, the statistical uncertainties for each model does not 
vary significantly if the overall coupling is allowed to vary 20-30%. 
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Table 1: Fermion couplings to the Z' for the considered models. An overall e/ cos^vk has 
been factored out. 
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3 Basic observables 



If a Z' is discovered at the LHC, the next step will be to determine the underlying model 
from which it arises. The following observables can be used to check whether the data fits 
the hypothesis of a certain model and begin to measure Z' properties. 

• Z 1 mass and total width, Mz> and Tz>- 

One should be able to find a peak at the LHC from an excess of dilepton events. The 
location of the resonance determines the mass of the Z 1 . The width is determined 
by fitting the resonance peak to the Breit-Wigner form 1/ [(M^ — Mf,) 2 + M|,r|,]. 
The width is sensitive to Z' couplings to all final states, and can probe invisible decay 
modes. We assume that the Z' has no invisible decays besides to neutrinos in our 
example models. Our analysis is nearly width-independent, and our results will not 
differ provided the invisible width is not large enough to make the Z' too broad. 

• Cross section to e + e~, a. 

In defining the Z' on-peak cross section, we follow [12] and keep events within ±317 
of the resonance peak. Using this instead of a fixed value allows for more consistency 
between models in isolating the Z' from the other neutral gauge bosons, and renders the 
cross section less sensitive to the width chosen. A nearly width-independent quantity 
is aT. The total cross section provides a good first separation between models and 
gives an indication of overall coupling strength and leptonic branching fraction. 

• Forward-backward asymmetry, A FB . 

A FB measures the relative difference of forward-scattered events and backward-scattered 
events: 



where F(y) = ^ d cos 9 dy f° os , B(y) = f^dcos dy f c a os9 , y is the Z' rapidity, and y max 
is the maximum allowed Z' rapidity given by |ln(s/M§,). The electron-quark angle 9 
is taken to be in the Collins-Soper frame [13], but there is an ambiguity in the quark 
vs. anti-quark direction, since it is unknown which proton carried it. We follow the 
suggestion in [14] and choose the quark direction along the direction of the Z' rapidity. 
Equivalently, one could choose one beam as the quark direction and and exploit the 
antisymmetry in y [15]. The value of y\ can be chosen to throw away events with low 
Z' rapidity, i.e., those where the quark direction is more likely to be misidentified. We 
study the dependence on y x in our analysis. A FB is quite sensitive to models with 
parity violating couplings. The on-peak value of A FB is defined by keeping events 
within ±3r of the resonance peak, as for the cross section. 
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Rapidity ratio, R. 

The central/forward rapidity ratio is defined as 



r: yi [F{y) + B{y)]dy 

V1 if;r+JZJ™ +B (yWy' 

R measures the ratio of central rapidities to extreme rapidities. Since the up and down 
PDF distributions have substantially different profiles, they should weight Z' events 
differently in rapidity, and thus R can help distinguish up versus down couplings. R 
is defined by keeping events within ±3r of the resonance peak. 

Off-peak asymmetry, A°J^~ peak . 

In addition to the above on-peak observables, the profile of Afb in dilepton invariant 
mass bins below Mz< may vary considerably among models. For this observable, we 
integrate instead in the region 2/3Mz> < Mu < Mz< — 3IV. 



4 Calculation 

To study the Z' signal at the LHC, we perform a fully differential next-to-leading-order QCD 
calculation of all observables considered. We include all spin correlations and interferences 
between the photon, Z, and Z' . We impose the following basic acceptance cuts on the final 
state lepton transverse momenta and pseudorapidities: p l T > 20 GeV and \r) \ < 2.5. Previous 
studies have found that detector resolution effects and other measurement errors are unlikely 
to have a significant effect on the e + e~ final state [12], and are neglected. A CMS simulation 
of Z' production found reconstruction efficiencies near 90% in the electron channel and no 
significant detector systematic errors [16] . In addition, electron energies can be measured to 
better than 1% accuracy, and invariant masses can therefore be reconstructed very well. The 
factorization and renormalization scales are taken to be fip = /ir = Mz>- They are varied 
simultaneously from M^/2 to 2Mz> to determine scale errors. We use the CTEQ 6.5 NLO 
PDF sets [17]. PDF error estimates are determined by calculating each observable with each 
of the 40 PDF eigenvector sets (corresponding to 20 +/- directions in parameter space), and 
combining the errors in each +/- shift in quadrature. Statistical errors are those that can 
be expected for an integrated luminosity of 100 ftT 1 unless stated otherwise. 

We present below in Tables[2j[3j and|4]for orientation the cross section, cross section times 
width, acceptance, asymmetries, and central-forward ratio for our four example models. We 
examine Z' states with Mz> = 1.5 TeV and 3 TeV. The acceptance denotes the fraction of 
events that pass the cuts on p l T and rf presented above; the fully inclusive results are obtained 
by dividing the cross section results by the acceptances. There are a few interesting features 
to note in these numbers. 

• The acceptances after imposing realistic cuts are independent of the model considered 
to the percent level. This was also observed in a simplified analysis in [10]. 



5 



• The acceptance, which is formed from a ratio of the cross section including cuts over the 
inclusive cross section, has tiny residual scale errors, indicating that it is insensitive to 
uncalculated higher order QCD corrections. We have checked that the leading-order 
acceptance is nearly indentical to the NLO result. It is generically true that scale 
errors are negligible for ratios of sufficiently inclusive quantities, and we see in Tables [3] 
and [4] that it is true for the other quantities considered here such as R, Apb, and 
A°Jg~ peak . This indicates that NNLO QCD corrections have no effect on the analysis 
of Z' properties. 

• While yi = 0.8 is the canonical choice for Afb measurements, this is an extreme 
rapidity for a 3 TeV Z' at the LHC Therefore, we also give 3 TeV values for yi = 0.4 
in Table H 

• PDF errors are relatively large for the total cross section, and are not negligible for 
other observables. 

• NLO results are substantially different from LO, especially for the total cross section. 
For instance, at LO, for the \ model at 1.5 TeV, one finds a = 40.0 fb ±23 and 
A°/ B = -0.2137±°;°°°| (errors are scale) versus 50.3 fb and -0.217 at NLO.' NLO 
corrections are larger than LO scale errors would suggest. This is not surprising; 
similar results have been seen for Z production for the scale range considered here [18] . 
NNLO corrections to the Z cross section have been shown to leave the NLO central 
value essentially unchanged while further reducing the scale error, and so we have 
confidence in our NLO analysis. 

Before continuing, we note a caveat regarding our calculation. We have not included 
higher-order electroweak effects. The complete 0(a) electroweak corrections for SM 
production for hadron colliders have been calculated in [19]. Several components of the 
higher-order electroweak corrections can be identified, and we indicate below how we expect 
them to affect our study. 

• QED corrections can be separated into two classes, those associated with initial and 
final state radiation. Those associated with initial state radiation lead to collinear 
singularities that must be absorbed into the bare PDFs, and modify the DGLAP 
evolution of the PDFs. The initial-state QED corrections in a DIS scheme were shown 
to be at the percent level or smaller [19]. This study did not include the QED effects on 
PDF evolution, as the appropriate PDF sets were unavailable at the time. The initial- 
state QED effects on DGLAP evolution have since been incorporated into a global fit 
to the available data [20], and they appear unlikely to have a significant affect on our 
analysis. 

• Final state QED radiation can have a significant effect on the lepton pair invariant 
mass distribution [19,21]. They should be included in a more complete analysis. 
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• Weak corrections consist of factorizable terms that can be absorbed into effective cou- 
plings and masses, and non-factorizable pieces arising from box diagrams that cannot. 
Since we treat the Z 1 couplings to fermions as free parameters in our coupling extrac- 
tion, the factorizable corrections have no effect on our study. They enter only when 
we choose a value of sin 2 9 W for the left-right model results. 

• The non-factorizable corrections arising from electroweak logarithms of the form lnf^f — ) 

lvl w,z 

become large for high lepton pair invariant masses in the Standard Model. While we 
expect these contributions to be small on the Z' peak, they can become important 
in off-peak observables. The new box diagrams are those containing both a Z' and a 
Z. These should contain only a single logarithm [22], and we expect them to be less 
important than for Standard Model lepton pair production at high energies. 



5 Measuring charges 

Ideally, one would like to be able to measure the couplings directly, rather than refer to 
a particular model. We study the extraction of four coupling combinations that can be 
determined from the on-peak observables considered above. Two of the combinations are 
the parity symmetric combinations c u , Cd introduced in [10], which we introduce explicitly 
later. We extend that parametrization to study the extraction of parity violating coupling 
combinations. 

First, we note that the differential cross section can be expressed as 

= £ + + e i) + ^ ~ - <£)]■ ( 7 ) 

q=u,d 

We have assumed that the couplings are generation- independent, thus the contributions from 
different quarks of the same type can be combined. The coefficients a\ and a\ represent 
what is left after the couplings are extracted from the differential cross section. They are 
composed of the PDFs and matrix elements integrated over phase space subject to the cuts 
discussed above. We have checked that the 7 and Z interference and squared terms are 
negligible, and we have dropped them in this parameterization. 

Let us define, as in [10], the parity symmetric coupling combinations 

c q =^(<1 2 R + dM + 4) = K + d)Br(Z' - e+e-). (8) 

We further define the parity-violating combinations 

M z , 
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Table 2: Cross sections and acceptances for our example Z' models. From top to bottom, 
the errors are statistical (100 fb -1 ), PDF, and scale. Errors in a ■ T are those due to the cross 
section only. As the acceptance is a theoretical quantity, we do not report a statistical error. 

which can be accessed via measurements of the differential cross section. Our equation for 
the differential cross section now reads 

E -Yj-M^ + al^}- ( 10 ) 



dyd cos 9 Mz> 

q=u,d 



We absorb the overall factor into the coefficients, a\ 2 = j!^a\ 2 , so that 

djj^e = Z> [alc * + a * e * ] - 

q=u,d 

In the narrow width approximation, the a\ 2 scale as a\ 2 ~ 1/r. Switching to a\ 2 removes 
almost all width dependence from these factors; we have checked that they vary by less than 
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Table 3: Forward-backward asymmetries and central to forward ratio for our example Z' 
models. From top to bottom, the errors are statistical (lOOfb -1 ), PDF, and scale. 

0.5% over the range of widths considered here. The only dependence of the a\ 2 on the Z' 
being considered is through Mz>- These coefficients therefore need be computed only once 
for a given Mz> and set of cuts. We present them below for Mz> = 1.5 TeV. 



We propose to use Eq. 11 to determine the four quantities c q and e q . By integrating Eq. 
11 over four different regions in y and 9, one obtains four equations for the four unknowns. 
The left-hand-side is determined by experiment and the integrated coefficients a±, a<i are 
determined theoretically. Thus one can solve for the unknown c q and e q . While in principle 
any four independent regions could be used, we would like to minimize the errors by isolating 
each of the four unknowns as much as possible in our four equations. To do this, we use 
the four observables F< = J\ dy F (y) , B < = J!^dyB{y),F > = (g ma * + ^ ym JdyF{y), 
and B > = ( f Vmax -f f Vl )dyB(y). Comparing F and B helps to separate c from e : one 

J Vl J Umax' ' " ' ' y y 

expects c q to contribute to F + B and e q to F — B. Separating different Z' rapidities helps 
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±° 


0.152 


_i_0.007 
=^0.009 

±° 


1.138 


i 0.068 
=1 -0.076 
_i_ 0.006 
=1= 0.004 



Table 4: Forward-backward asymmetries and central to forward ratio for yi = 0.4 and 
M z > = 3 TeV. 



to isolate up-type from down-type couplings due to their different PDFs. These quantities 
are related to the forward-backward asymmetry, rapidity ratio, and cross section via 



2 U ApB R yi + l h 



2 R yi + 1 



2 V R 



.VI 



1 



(12) 
(13) 
(14) 
(15) 



After integrating Eq. 11 over these four regions, we end up with the system of linear 
equations 



rh = Mc 



(16) 



where 



m 



( F < \ 



( c u \ 



(17) 



\ B > J 

and M is a matrix composed of the coefficients a\ and a\ integrated over the appropriate 
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ranges of y and 9. Explicitly, M takes the form 





a\ 


S,< 


af 


U< 


«2 


u< 


<4\ 


Ib < 


a L l 


Ib < 


af 


Ib < 


a u 2 


Ib < 


4 


Jf > 


a\ 


If > 


af 


If > 


a\ 


If> 


aj 


V s B> 


a'[ 


Jb > 


af 


Ib > 


«2 


Jb > 


4 ) 



M 



The entries are determined by running our code for certain basis models with charges such 
that one of the c q or e q are equal to 1, and the others are zero. The photon and Z contributions 
have been checked to be negligible and are turned off in this process. For example, with 
Hi = 0.8, using the central PDF, we get 



M 



/ 5638 


4175 


1747 


828 


\ 


5638 


4175 


-1746 


-827 




3610 


1519 


2101 


784 




^ 3610 


1519 


-2101 


-784 


/ 



fb 



(19) 



for Mz> = 1.5 TeV. Solving for the couplings is now straightforward: c = M^m. As 
mentioned we have found that these entries vary by less than 0.5% over the range of widths 
for the examined models. 

Below we illustrate this extraction procedure applied to our four example models. We 
simulate experimental results for F <1 B <) F>, and By using the full NLO cross section 
including the 7 and Z contributions. We assume an integrated luminosity of 100 fb -1 . We 
then subject these results to our extraction procedure described above. Doing so, we obtain 
"measurements" of c q and e q . These measurements are compared to the theoretical input 
parameters used in generating the data. Our results are shown below in Table [5j For 
Mz' = 1.5 TeV we use y\ = 0.8; for 3 TeV we use y\ = 0.4. We have checked that scale 
errors are negligible in this analysis, as discussed previously, and they have been omitted. 

We see that there is good agreement between the couplings extracted using this method 
and the theoretical input values. One caveat must be discussed. Clearly, our simulation of 
the experimental results is simplistic. However, it allows us to study two important issues. 
First, it shows us that this technique is consistent, and particularly that the neglect of the 7 



and Z terms in Eq. 11 is justified. It also allows us to study how accurately the c q and e q can 
be determined at the LHC given the expected theoretical and statistical errors. For 100 fb -1 , 
individual errors on c q and e q can be rather large, especially in the 3 TeV case. However, we 
will see that when the couplings are taken together, models are well-discriminated, especially 
at the SLHC. 

We now illustrate these results graphically for Mz> = 1.5 TeV, first for the LHC assuming 
100 fb" 1 and then for the SLHC assuming lab . The LHC results for c uc i and e u d are 
projected into a c^-plane and e^-plane, and shown in Figs. [T] and [2j respectively; SLHC 
results are in Figs. [3]and|4j We also present results for Mz> = 3 TeV for the SLHC in Figs. [5] 
and |6j the statistical errors at the LHC are simply too large for a meaningful coupling 
extraction without more data. In general, the errors are ellipsoids in the four-dimensional 
space of (c u , cj, e u , e^). The statistical errors are very nearly diagonal in c u /cd vs. e u /ed due 
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M z < 


Mdl 


c u x 10 3 


c d x 10 3 


e„ x 10 3 


e d x 10 3 






n 66 


3 30 


n 


—2.11 




X 


O 68 ± 14±2oo 

U.UU — 1 — W • X t: — 1 — r\ oo 


3 32 ± 24±n^ 


01 ± 76±Hi 

u.ui _i_ u. i u_i_q 35 


— 2 09 ± 1 92±?™ 

jLj • \J tj _l — X . • V — — 1 — 1 r 7C\ 




rl> 


0.81 


0.81 








1.5 


0.81 ± 0.10=bg;g| 


0.81 ± 0.17=bg;J§ 


0.01 ± 0.53±[j- 01 


0.00 ± 1.32±g;gl 


TeV 


V 


1.08 


0.67 





-0.24 




1.09±0.11±Hq 

u.uy 


0.68±0.18±S-?2 

U. 10 


0.01 ± 0.57±S-85 


-0.23±1.41±8-92 

u. iy 




LR 


1.59 


2.69 


0.60 


1.03 




1.61 =b 0.16d=2-lg 


2.71 ± 0.26±^i? 


0.57 ± 0.84±2-$ 


1.12 d= 2.10d=8-f| 




X 


67 


3 36 


n 


—2.15 




69 ± 1 4d=9-?i 

U .U>7 _l — X • Jl — 1 — 1 1 C 


3 36 ± 3 02±? oo 

"7 . iai _i — UiUai _ I — i no 


00 ± 3 26±S o? 


— 2 12 ± 7 81±?'?n 

— . X — 1 1 . K..' X 1 1 ■ > 




1> 


0.82 


0.82 


o 


o 


3 


82 ± 1 01±Ho 


82 ± 2 22±R-5n 


01 ± 2 36±n 


—0 00 ± 5 62±S 


TeV 




1.10 


0.69 





-0.25 




V 


1.09±1.10±°;|| 


0.69 ± 2.41±°|i 


0.00 ± 2.55=bg-g| 


—0.25 =b 6.07=bg;S| 




LR 


1.62 


2.74 


0.61 


1.05 




1.63±1.6±°;j» 


2.74±3.45±J;$ 


0.57±3.68±g;J| 


1.09±8.78±°;| 



Table 5: Result of extracting Z' couplings at the LHC. Theoretical values for the couplings 
are listed first. The results of our extraction procedure are shown next. Errors are statistical 
and PDF, respectively. 

to our choice of measurements. However, the PDF and therefore the combined errors are 
not. For simplicity of presentation we plot the projections of the 4-D PDF and combined 
errors into the c Ujd and e U)d planes. We also conservatively take the larger of the +/- PDF 
deviations for each ellipse. We have plotted a contour for the family of Eq models on the 
M' z = 1.5 TeV plots for reference. We note the following trends in these plots. 

• The statistical error ellipses are rather narrow, with minimal extent in the c u + q and 
e u + ed directions. This occurs because statistical errors are reduced when we add 
contributions form up and down quarks, rather than attempt to distinguish between 
them. 

• The PDF error projections are nearly orthogonal to the statistical errors, and are 
minimized in the c u — q and e u — directions. We see that models are easily distin- 
guishable with only statistical errors due to the narrow ellipses, but when including 
PDF errors, a large volume of coupling space is occupied. However, for a 1.5 TeV Z', 
the c couplings can still distinguish models. One can still determine these coupling 
combinations within a reasonable window. 
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• The PDF errors scale with the couplings. For models with zero values, such as the ijj 
model for e u and e^, the errors are mostly statistical. 

• The e couplings are harder to measure. Much of this stems from the fact that the 
difference F — B is used to extract them, which has higher statistical error than the 
corresponding F + B for the c couplings. 

• All Eq models have e u = and < 0. A substantial departure from this could rule 
out this family, though due to large statistical errors the SLHC might be needed. 

• The error ellipses are quite narrow in Figs. [3] and |4| for a 1.5 TeV Z' at the SLHC, 
allowing reasonable determination of the couplings. The errors will improve further 
should PDF errors improve. 

• Charge extraction for M Z i = 3 TeV is difficult, even at the SLHC with present PDF 
error estimates. 



Z' c Couplings, y=0.8, M z ,= 1.5 TeV 

0.0020 



0.0015 



? 0.0010 



0.0005 



0.0000 

0.000 0.001 0.002 0.003 0.004 

c d 

Figure 1: Simulated measurements of the c u ^ couplings at the LHC for our test models. The 
dashed ellipses are the statistical errors expected for M z > = 1.5 TeV and lOOfb -1 of data, 
the dotted ellipses are the current estimated PDF errors, and the solid ellipses denote the 
combined errors. The Eq family of models lie on the dot-dashed contour. 

5.1 Distinguishing Models 

We now test whether the above observables are sufficient to distinguish our example models 
at the LHC. We compare models pairwise by assuming that one is correct and finding the 
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Z' e Couplings, y=0.8, M z ,= 1.5 TeV 



0.002 



-0.001 









— 




— 


— \ 7 '^* 













-0.004 -0.003 0.000 0.002 0.004 

e d 

Figure 2: Simulated measurements of the e u ^ couplings at the LHC for our test models. The 
parameters and ellipses are as discussed in the previous plot caption; the E 6 contour is a 
line segment. 



Z' c Couplings, y=0.8, 1 ab \ M z ,= 1.5 TeV 

0.0020 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




0.0000 1 1 1 1 ' 1 1 ' ' 1 1 ' 1 1 1 1 1 — "t- 1 ■ f 1 

0.000 0.001 0.002 0.003 0.004 

C d 

Figure 3: Simulated measurements of the c u ^ couplings at the SLHC for our test models. 
The dashed ellipses are the statistical errors expected for Mz> = 1.5 TeV and 1 atT 1 of data, 
the dotted ellipses are the current estimated PDF errors, and the solid ellipses denote the 
combined errors. The E$ family of models lie on the dot-dashed contour. 
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Z' e Couplings, y=0.8, 1 ab 1 , M z ,= 1.5 TeV 

0.002 



0.001 



J 3 0.000 



-0.001 



-0.002 

-0.004 -0.002 0.000 0.002 0.004 

Figure 4: Simulated measurements of the e u ^ couplings at the SLHC for our test models. 
The parameters and ellipses are as discussed in the previous plot caption; the E 6 contour is 
a line segment. 



Z' c Couplings, y=0.4, 1 ab \ M z , = 3 TeV 

0.003 



0.002 



? 0.001 



0.000 



-0.001 

0.000 0.002 0.004 0.006 

C d 

Figure 5: Simulated measurements of the c u ^ couplings at the LHC for our test models. 
The dashed ellipses are the statistical errors expected for Mz> = 3 TeV and 1 ab" 1 of data, 
the dotted ellipses are the current estimated PDF errors, and the solid ellipses denote the 
combined errors. 
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Z' e Couplings, y=0.4, 1 ab \ M z , = 3 TeV 

0.002 p 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 




-0.002 u ' 1 1 1 1 ' 1 1 ' 1 ' 1 1 ' 1 ' 1 1 ' 1 ' 1 1 u 

-0.004 -0.002 0.000 0.002 0.004 

Figure 6: Simulated measurements of the e u ^ couplings at the SLHC for our test models at 
Mz> = 3 TeV. The parameters and ellipses are as discussed in the previous plot caption. 

chi-squared of the other model as a test, using the errors of the first model. Statistical and 
PDF errors are combined in quadrature. While the PDF errors roughly represent a 90% 
confidence level, the precise meaning isn't clear [23]. We conservatively take them to be la 
in the combination and again take the larger of the +/- deviations, as in our plots. Scale 
errors are unimportant in this analysis and have been dropped. 

We use the charges c q and e q , as well as A°Jl~ peak ' 0A to form the x 2 f° r 1-5 TeV. Rather 
than use the charges c q and e q directly, we form the x 2 by diagonalizing the errors in these 
variables, as in our plots. This allows us to best exploit the model separation in directions 
where the errors are minimized. Afb numbers tend to contribute the most to x 2 since their 
errors are relatively small, and this is the only off-peak data we have included here. If one 
uses only ApBi choosing a lower rapidity gives more data, reducing the error. However, for 
too low a rapidity, quark direction misidentification tends to wash out the actual values. We 
have found that y\ = 0.4 is a good choice when not including other data. 

We see from Table [6] that for Mz> = 1.5 TeV, these observables should distinguish the 
considered models quite reliably. As values of 1 — 2 a correspond to confidence levels of 68% 
and 95% respectively, we see that most are distinguishable with 99% C.L. or greater. While 
the ip an d Tj models have similar values for all observables, the errors are small enough that 
together one can separate them with a confidence level of nearly 90%. The other models are 
easily distinguishable with 100 ftT 1 of data. In the full c/e analysis, yi = 0.8 appears to be 
the most discriminating choice. For the heavy 3 TeV Z', there are very few high rapidity 
events; we have therefore restricted our analysis to y\ = 0.4. However, there is still some 
distinguishability at 100 ftT 1 with this choice, despite the large errors seen in the plots and 
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Table 6: Pairwise y 2 values for our model comparison, for 100 fb . The separate y 2 con- 
tributions from the on-peak c g , e q couplings have been shown, as have the total y 2 values 
including A°ff~P ea > 0A ? anc j the corresponding confidence that the models are distinct in stan- 
dard deviations. The model in each row is assumed to be the correct, measured model, and 
is tested against the hypothesis in each column. We have performed this test for several 
choices of on-peak y x ; for M z > = 1.5 TeV, y\ = 0.8 appears to be the optimal choice. For 
the off-peak asymmetry, only y\ = 0.4 has been used. For Mz> = 3 TeV only y\ = 0.4 has 
been used. Note that since the statistical errors come from the row models, and PDF from 
the columns, this table is not symmetric. 



in Table [5j Several models can be separated at 1 — 2a (68-95% CL). We see below in Table [7] 
that the situation improves signficantly with the SLHC, with all models but the ip/i] pair 
well-separated. 

In principle, one could use the off-peak data in a full analysis of the charges by modi- 
fying Eq. [7] and adding the neglected Z and photon contributions. The expression for the 
differential cross section becomes 



d 2 a 



dyd cos 9 



Yl + 4)(4 + 4) + 4(Qr ~ ?D(4 - 4) + b\q R e R + b q 2 q R e L ] 



=u,d 



+ hq^R + biq L e L + c. (20) 

The known Z and photon charges have been folded into the coefficients b and c. There are 
five unknowns, and R , but only four independent combinations of these appear above. 
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Xle Xtot a 
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37 43 5.5 
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15 29 4.3 




1.1 2.3 0.2 


4.6 26 3.9 


rj 


15 22 3.4 


1.3 2.3 0.2 




13 24 3.7 


LR 


14 14 2.4 


44 58 6.7 


30 38 5.1 





Table 7: Pairwise x 2 f° r 1 & b \ Hi — 0-4, and M z > = 3 TeV. As before, the rows are tested 
against the hypothesis columns. 

This can be seen by noting that increasing the quark charges by a factor of two and decreas- 
ing the lepton couplings by a factor of two leaves the differential cross section unchanged, 
indicating that a degeneracy exists. Modifying our above procedure would require solving 
four quartic equations for the four independent unknowns. However, the linear terms in q x e 
provide information on the signs of the charges, which is an enticing prospect; it is likely one 
can gain significant statistical precision, and determine the signs with confidence. In addi- 
tion, if one can determine the invisible width, this would also yield Br(Z' — > e + e~), which 
could be used to separate q x e and solve for the individual charges qi,R and cl,r-, with signs 
if off-peak data is analyzed. The ILC should also be able to probe the interference region 
for the masses considered [24], and may also be able to determine directly, breaking our 
q x e degeneracy. 

6 Conclusions 

In this paper we studied the measurement of Z' couplings at the LHC We performed a fully 
differential NLO QCD calculation of the Z' signal with all spin correlations and interference 
effects with the SM 7 and Z included, as well as realistic LHC acceptance cuts. Using four 
example models arising from grand unified theories, we quantified the effect of statistical, 
PDF, and residual scale errors on important Z' observables. We found that residual scale 
errors are negligible for observables formed from cross section ratios, such as Aps, indicating 
that corrections from NNLO QCD effects are unimportant in Z' studies. However, statistical 
and PDF errors are significant and have approximately equal effects with lOOfb -1 at the 
LHC. 

We introduced a set of Z 1 coupling combinations that can be determined from on-peak 
measurements at the LHC. The parity symmetric combinations c u ^ which can be accessed by 
measurements of the inclusive cross section were introduced previously in [10]. We extended 
this parametrization to include the parity violating combinarions e Ut d that can be probed 
once differential measurements are made. The differential cross section factors into a sum 
over products of these couplings times transfer functions that depend on the model under 
consideration only through the Z 1 mass; the slight dependence of these functions on the 
width was found to be less than 1%. These transfer functions only need to be evaluated 
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once for a given Z' mass and set of cuts, and then can then be used in simulations regardless 
of the underlying Z' model. They form a matrix which connects measurements in different 
kinematic regions to the underlying Z' couplings. To access these transfer function one 
needs only to run a Z' simulation code for basis vectors in coupling space of the form 
(c u , c d , e u , e d ) = (1, 0, 0, 0), etc. 

We computed these transfer functions and used them to examine how well Z' couplings 
can be determined at the LHC assuming 100 ftT 1 , and at the SLHC assuming lab -1 . As 
illustrative examples we again used the four example models discussed previously. Both 
statistical and PDF errors give equally important contributions to the uncertainty in cou- 
pling measurements at the LHC. The statistical and PDF errors turn out to be maximal in 
orthogonal directions in both the (c u , q) and (e u , e^) planes. We found that the c u ^ could 
be determined with reasonably good precision at the LHC, easily well enough to distinguish 
between the four example models considered. However, the e u ^ will be relatively poorly 
determined with 100 fb" 1 assuming current PDF errors; a more accurate measurement of 
these couplings sufficient to tell which of the four example models they came from required 
1 ab _1 in our analysis. 

Our analysis can be extended in several ways. Inclusion of off-peak observables can 
improve the precision of the coupling extractions; we saw in Section [6] that the off-peak 
asymmetry increased the \ 2 i n the point-wise comparison between models. However, off the 
Z' peak interference terms between the Z' and SM 7 and Z give the dominant contributions. 
Our coupling parametrization must be enlarged to include these effects. Another direction 
in which to extend our analysis would be to break the degeneracy between lepton and quark 
couplings to the Z' . This degeneracy is clearly visible in Eqs. [7] and |20j if we scale the quark 
couplings up by a factor of two; scaling the lepton couplings down by the same factor leads to 
an unchanged cross section. Measurements of rare Z 1 decays can break this degeneracy [25] . 
In principle the Z 1 width could break this degeneracy if the invisible decay width was known. 
The Z 1 invisible width is also of interest in models where the Z 1 acts as the messenger to a 
hidden sector, as it gives some insight into the matter content of this sector of the theory. 
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